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Abstract While clusters of galaxies are considered one of the most important cosmo- 
logical probes, the standard spherical modelling of the dark matter and the intracluster 
medium is only a rough approximation. Indeed, it is well established both theoretically 
and observationally that galaxy clusters are much better approximated as triaxial ob- 
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jects. However, investigating the asphericity of galaxy clusters is still in its infancy. We 
review here this topic which is currently gathering a growing interest from the cluster 
community. We begin by introducing the triaxial geometry. Then we discuss the topic 
of deprojection and demonstrate the need for combining different probes of the cluster's 
potential. We discuss the different works that have been addressing these issues. We 
present a general parametric framework intended to simultaneously fit complementary 
data sets (X-ray, Sunyaev Zel'dovich and lensing data). We discuss in details the case of 
Abell 1689 to show how different models/data sets lead to different haloe parameters. 
We present the results obtained from fitting a 3D NFW model to X-ray, SZ, and lens- 
ing data for 4 strong lensing clusters. We argue that a triaxial model generally allows 
to lower the inferred value of the concentration parameter compared to a spherical 
analysis. This may alleviate tensions regarding, e.g. the over-concentration problem. 
However, we stress that predictions from numerical simulations rely on a spherical anal- 
ysis of triaxial halos. Given that triaxial analysis will have a growing importance in 
the observational side, we advocate the need for simulations to be analysed in the very 
same way, allowing reliable and meaningful comparisons. Besides, methods intended to 
derive the three dimensional shape of galaxy clusters should be extensively tested on 
simulated multi-wavelength observations. 

Keywords Cosmology • Galaxy Clusters • Triaxiality ■ Gravitational Lensing ■ 
X-rays • Sunyaev Zel'dovich 



1 Three-Dimensional Shape of Galaxy Clusters 

Spectroscopic galaxy redshift surveys and numerical N-body simulations have revealed 
a large scale distribution of matter in the Universe featuring a complex network of 
interconnected filamentary galaxy associations. Vertices, i.e. intersections among the 
filaments, correspond to the very dense compact nodes within this cosmic web where 
one can find massive galaxy clusters. 

In this review, we concentrate on the shape of galaxy clusters. 



1.1 Galaxy Clusters are not Spherical 



There is much observational evidence for clusters not being spherical objects, from 
the non-c ircular projection of variou s probes: optica l, from the density maps of cluster 
galaxies I Carter and Metcalfe 1980l: Binggeh IQSd): X-ray, f r om the surface bright- 



ness maps ( Fabricant et al.l 1984 ; Buote and Canizares 199^ 19961: Kawaharal I2OI0I : 



iLau et al]|2012l ); SunyaevZel'dovich pressure maps (|Savers et al.l 2011al'): strong gravi- 



tational lensing dSoucail et al.|[l987l ). and weak gravitational lensing ()Evans and BridQ 
I2OO9I : lOguri" et al. I I2OIOI . I2OI2I ). Recently, the azimuthal variation of galaxy kinematics 
has been detected for the first time in a stacked sample of 1 743 galaxy clusters from 
the SDSS dSkielboe et al.ll2012h . They find that the line of sight velocity dispersion of 
galaxies lying along the major axis of the central galaxy is larger than those that lie 
along the minor axis. This detection provides further evidence for the asphericity of 
galaxy clusters. 

On the numerical side, haloes forming in cosmological simulations hav e been found 



to be triaxial in shape, with a preference for prolateness over oblateness (|F Venk_et_al 



19881 : iDubinski and Carlberg||l99ll : IWarren et a"l]|l992l : ICole and Lacevll996l : ljing and Sutj 
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2002':'Hop kins et al.l2 005":'Bai lin and Steinmetzl2005l:lKasun and Evrardl2005l:IPaz et al.l 
a)06a; Al kood et al.i ,200a : iBett et al.1 l2007l : iMunoz-Cuartas et all I2OI1I : ICao et alF 
2012). These simulations also predict an evolution of the shape with mass and red- 
shift: low mass haloes appear more sphe rical than high mass haloes (see, however, an 
opposite conclusion bv lRossi et al.ll201lh . and for a given mass, lower redshift haloes 
are more spherical than high redshift haloes. 

However, more than forty years after the first observational evidence of the as- 
phericity of galaxy clusters, the majority of cluster studies use the spherical assump- 
tion. Historically, this was due to the fact that the quality of the data may not allow a 
triaxial model to be constrained. Besides, studies often rely on the analyses of a single 
data set {e.g. lensing. X-ray, SZ or dynamics of cluster members only), while to recover 
triaxiality one needs to combine data sets, unless some priors are used (see Section 2). 

Besides the expectation for triaxial collapse coming from first principles (see Sec- 
tion [L2]), another hint regarding the need for non-spherical models is the mass discrep- 
ancy found between, e.g. lensing and X-ray data when spherical symmetry is a ssumed 
i Piffaretti et"ai]|2003l ; [clowe et al.|[20o3 : ICavazzilboOSl : ICorless and Kinel[2007[ ). 



1.2 Triaxiality: a Consequence of Gravitational Collapse for an Initial Gaussian 
Random Field of Density Fluctuations 

Interestingly, triaxial collapse is a straightforward prediction of structure growth driven 
by self-gravity of Gaussian density fluctuations. Well into the mildly non-linear regime, 
gravitational structure formation can be described by the Zel'dovich approximation, 
which models the motion of dark-matter particles as inertial motion in a suitably 
adapted time coordinate. This assumption defines a map between initial, Lagrangian, 
and final, Eulerian coordinates of any particle. The Jacobian of this map is called 
the Zel'dovich deformation tensor, here abbreviated as F. The matter density is then 
simply given by the inverse determinant of F, times the initial density. 

Under the (reasonable) assumption of an initially irrotational flow, a velocity po- 
tential ip exists in Lagrangian space whose gradient is the Lagrangian velocity field. 
The Zel'dovich deformation tensor can be represented by the unit matrix plus the 
Hessian (curvature) matrix of the velocity potential, 

Fij = Sij + didjip . (1) 

It is thus symmetric and has three real eigenvalues (1 -I- Aj), i = 1,2, 3. 

Since determinants are invariant under orthogonal transforms, the inverse determi- 
nant of the Zel'dovich tensor is 



(det FY 



(2) 



Collapse will set in according to the Zel'dovich approximation whenever any of the 
eigenvalues approaches —1. If all eigenvalues A^ were the same, collapse was spherical; 
if they were all different, col lapse was triaxial. 

It was shown already bv IPoroshkevichl ( 1970l ) that the probability distribution for 



the Xi in a Gaussian random field is proportional to the product of absolute differences 

|A,-A,| (3) 
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for any pair i,j with j ^ i. Tiiis means that equal eigenvalues are not allowed in a 
Gaussian random field, thus excluding spherical and even spheroidal collapse: in any 
realistic triple of eigenvalues, no pair of them can be equal. Collapse will proceed first 
along the principal axis belonging to the largest eigenvalue, forming the sheet-like 
structures called 'pancakes' by Zel'dovich. Those will then contract along the principal 
axis of the second-largest eigenvalue, forming essentially one-dimensional bridges of 
matter that finally shrink to triaxial haloes along their remaining axis. This inherent 
and necessary triaxiality of the mildly non-linear collapse is retained through the non- 
linear collapse of the haloes which is stopped by virialization. 



1.3 Triaxiality; an Outgrowth of the Large Scale Structure Formation Scenario 

The standard cosmological framework (ylCDM), which consists of a cosmological con- 
stant and cold dark matter (DM) with Gaussian initial conditions, envisages structure 
formation as a hierarchical merging process. With this perspective, gravity is con- 
stantly pulling lumps of matter together to form increasingly larger structures. The 
structures we see in the Universe today (galaxies, clusters, filaments, sheets and voids) 
are predicted to have formed in this way, with galaxy clusters sitting atop this hierar- 
chy and being the largest virialized structures that have had time to collapse under the 
influence of their own gravity. In particular, being a tracer of the primordial density 
perturbations, this large scale structure scenario leads to a picture where the matter 
is distributed as a network of gigantic dense (filaments) and empty (voids) regions, 
creating a vast foam-like structure called the "cosmic web", with the densest regions 
of the dark matter cosmic web hosting massive clusters of galaxies. Numerical sim- 
ulations, in agreement with analytic predictions, indicate that the infall of material 
into the most massive dark matter haloes is not spherical but is expected to be pref- 
erentially funnelled through the filaments where the haloes are embedded. The cluster 
mass haloes would indeed acquire most of their mass from major mergers along the 
filaments, hence leading to an alignment between the major axis of the host haloe 
and the large-scale filament (iBai lin and Steinme t 2JI2OO5I : lAltav et al.ll200^ : IPatiri et al.l 



I2OO6I : lAraedii -Calvo et al ] |2007l : lBrunino et al.ll2007l ). Clusters would then relax from 



this chaotic initial state to a quasi-equilibrium via violent relaxation. This process leads 
to equilibrium which is plausibly related to bounded triaxial DM distribution with a 
'universal' density structure, regardless of haloe mass, of cosmological parameters, and 
of the initial fluctuation spectrum. Therefore, triaxiality appears to be a direct out- 
growth of the large scale structure formation scenario, providing a record of the initial 
conditions in the Universe and of the topology of the cosmic structures. 



2 From 2D Observables to a 3D Mass Model: Deprojection 

2.1 An Under- Constrained Problem 

The problem of deprojection is typically under-constrained: we have access to 2D pro- 
jected informations (X-ray surface brightness maps, or lensing mass density maps) and 
aim to derive the 3D properties of the triaxial structure. 

As an example, let's consider gravitational lensing observations. Fitting triaxial 
models with lensing data only is an intrinsically under-constrained problem, since lens- 
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ing can never give full information about the 3D triaxial structure: it is constraining 
only the 2D projected mass density. This translates into much larger error bars on 
the parameters of a triaxial haloe constrained using lensing observations alone. These 
error bars fairly represent the true extent of our limited knowledge of the structure 
of galaxy cluster lenses, and make clear the importance of combining constraints from 
other theoretical or other observational data in order to narrow down the parameter 
space of a more realistic triaxial model. 

The first attempts to determine three dimensional morphologies were based on sta- 
tistical a££roaches consisting in the inversion of the distribution of apparent shapes. 
iHubblel (|l926l l first determined the relative frequencies with which galaxies of a given in- 
trinsic ellipticity, oriented at random, are observed as having various apparent projected 

ellipticities. Several following studies have then applied similar methods to different 

classes of astro nomical objec ts llNoerdlingeill979l:lBinggelilll98(]l:l iBinnev and de VaucouleursI 



1981; Fasano and Vio 



1991 



de Theiie et al 



1995; iMohr et alJ 



1 19951; 



Basilakos et al.l 



200d:ICoora v 2000':' Thakur and Chakrabortvi200ll ; jAlam and Rvdenl2002l ; lRvdenll99a 



Plionis et al...2004 ; et al.ll2006bir With the exception of disc galaxies, prolate-like 



shapes appear to dominate all cosmic structure on a large scale. 



2.2 Combining Data Sets 



We will consider three different types of data sets: gravitational lensing (both weak 
and strong), X-ray emission, and the thermal SZ effect. For details on these different 
data sets, we refer the reader to the other reviews published in this volume. 

When combining X-ray and SZ data sets, the idea is that one can infer the 3D 
properties of a cluster by taking advantage of the different dependences of the X-ray 
and SZ signals on the gas density and temperature: the SZ effect is proportional to the 
electron pressure integrated along the line of sight, whereas the X-ray surface brightness 
is proportional to the integral along the line of sight of the square of the electron 
density. Besides, gravitational lensing provides a direct probe of the two dimensional 
mass distribution projected along the line of sight. 

Combining complementary data sets to reconstruct the three dimensional proper- 
ties of galaxy clusters is not a new idea, and different authors have proposed different 
approaches. 



2.2.1 Theoretical Studies 



We review here the works that have developed methods which have been tested on 
simulated data sets onH 



IZaroubi et al.l (|1998l ) proposed a non parametric deprojection method specifically 
designed for the deprojection of X-ray, SZ and lensing maps of galaxy clusters, under 
the assumption of axial symmetry of the cluster. This method was first applied to a 
simple analytic model for cluster dark matter and gas distributio ns, then it was teste d 
on cosmological hydrodynamical simulations of galaxy clusters I Zaroubi et aLllioOll ). 
The authors found a good agreement between the actual (simulated) and reconstructed 
three-dimensional properties. 

Some studies proposed a non parametric deprojection method based on Abel's inte- 
grals. This inversion is appli ed on SZ and X- r ay da ta in order to infer temperature and 
density profiles of the ICM. ISilk and Whit3 l| 19781 ) were the first to propose to apply 
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the Abel inversion to the X-ray and SZ profiles, but their approach was limited to the 
estimation of the central values of gas density and temperature aiming at the deter- 
mination of cosmological parameters. More recentlv. lYoshikawa and Sutol 1 19991 ) used 
this technique for a non parametric reconstruction of radial density and temperature 
profiles using anal ytical and simulated cluster models. This method assumes spheri- 
cal symmetry, and lYoshikawa and Sutol ( 1999h mention the inclusion of non-spherical 
effe cts as an im portant next step. 

iReblinskvl !^00& l developed a non parametric algorithm for the simultaneous depro- 
jection of X-ray, weak lensing, and SZ data. They specify a geometrical model for the 
cluster assuming axial symmetry. Using gas-dynamical simulations, they demonstrated 
the qualit y of the deprojections. 

Later. IPuchwein and BartelmamJ 1 20061 ) proposed a method based on Richardson- 
Lucy deconvolution to reconstruct the three-dimensional gas density and temperature 
distributions in galaxy clusters from combined X-ray and SZ observations. They tested 
their algorithm against synthetic observations using both analytically and fully numeri- 
cally simulated clusters. They found that their method reconstructs the gas density and 
temperature distributions accurately in three-dimensions, even if observational noise 
is present. Additionally, they discussed a method to constrain the cluster inclination 
along the line of sight using X-ray temperature maps. They came to the conclusion 
tha t the method allows to reach a level of accuracy of ~ 15%. 

IFox and Pen! l|2002l ) considered the problem of deprojecting aspherical clusters. 
They first constructed a parametrised 3D axisymmetric cluster model, and determined 
the 3D cluster shapes using a fitting between the model predictions and the simu- 
lat ed data. 

iMarshaU et"al] (|2003l ) presented a Bayesian joint analysis of cluster weak lensing 
and SZ data, assuming spherical symmetry. This methodology was applied to two sets 
of simulated SZ and we ak lensing data sets. 

iLee and Sutol ( 20041 ) considered a deprojection method combining SZ and X-ray 
data and applied it to analytical cluster models, considering tnaxtal haloes with con- 
sta nt axis ratios. 

lAmeglio et al.l (|2007l . [20091 ) presented deprojection methods (parametric as well as 
non parametric) in order to recover the three dimensional density and temperature 
profiles by combining SZ and X-ray surface brightness maps, assuming the cluster to 
be spherically symmetric. They apply their techniques to a set of hydrodynamical sim- 
ulations of galaxy clusters in order to estimate the biases and scatters on the recovered 
masses^ 

Iau ison et al.l l|201lh proposed a parametrised spherical model of the intracluster 
medium aimed for jointly analysing SZ and X-ray data. This entropy based model which 
relies on the assumption of hydrostatic equilibrium is tested against mock observations 
of clusters from N - body/ hydrodynamic simulations. 

ISamsing et al.l ( 20121 ) presented a novel method for measuring a radially dependent 
shape along the line of sight of the ICM from the X-ray observations only. The method 
hinges on the assumption that the shapes, temperature and density profiles can be 
described by parametrised functions. This model generates fake spectra to be compared 
with the observed spectral information. A clear advantage of this approach is that it 
does not require any combination with independent measurements of e.g. the cluster 
mass or density profile. The major downside for the current observations is that it 
requires data of impressively high quality (~ 10® photons) to get a 5(t detection of the 
shape. 
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2.2.2 ObservaUonal Studies 

As seen before, from their different dependencies on the density, combining X-ray and 
SZ aUows us to directly infer the elongation of the gas distribution. No assumption is 
needed about hydrostatic equilibrium. The method exploits X-ray spectroscopic and 
photometric data plus measurements of the Sunyaev-Zel'dovich effect. One needs to 
deproject the X-ray and SZ data. Both data sets are supposed to trace the same 
temperature (remember that SZ gives a mass weighted temperature, unlike X-ray), 
and one assumes that no dumpiness and no contamination from structures along the 
line of sight might bias the SZ data. The gas distribution is modelled with an ellipsoidal 
parametric profile which can fit the observed X-ray surface-brightness and temperature. 
Comparison with the SZ amplitude fixes the e longation along the line of sight e^. For 



an isothermal plasma ( De Filippis et aDl2005h . 



where is the angular diameter distance to the cluster, SB is the surface brightness, 
ATgz the SZ temperature decrement, T the temperature, and Ax the emissivity. 
Finally, Bayesian inference allows us to deproject the measured elongation and the 
projected ellipticity measured in the plane of the sky, in order to constrain the intrinsic 

shape and orientation of the cluster. 

This approach was first employed in iDe Filippis et al.1 (|2005h . who considered a 



sample of 25 X-ray luminous clusters. They used parametric ellipsoidal profiles of con- 
stant eccentricities aligned along the line of sight, assuming an isothermal beta model 
parametrisation for the ICM. They found that the spherical hypothesis is strongly re- 
jected for most of the 25 clusters of the X -ray selected s ample studied. Considering the 
same data but assuming axial symmetry. et al.l ( 20061 ) showed that this sample 

is composed of a mixed population of prolate and oblate haloes, with prolate shapes 
preferred in ~ 60-76% of the cases. They observe an excess of clusters elongated along 
the line of sight, with respect to what is expected from a randomly oriented cluster 
population. They claim that a more general triaxial morphology might better describe 
the morphology of these clusters. Both studies acknowledged that adding gravitational 
lensing data would help further to constrain the shape. Indeed elongation strongly en- 
ters in the lensing properties of a haloe. Therefore gravitational lensing features brings 
inf ormation on the elongation of the total matter distribution. 

iMahdavi et al. I iioo?) provided a framework for the joint analysis of cluster obser- 
vations (JACO) which fits the mass models simultaneously to X-ray, SZ, and weak 
lensing data. They applied this method to Abell 478, assuming spherical symmetry. 
They do find a good agreement between all data sets, which points out that the spher- 
ical symmetry hypothesis was well motivated in the case of Abell 478. 

However, we stress that it can be questionable to combine different data sets assum- 
ing spherical symmetry; if a given galaxy cluster exhibits a significant mass discrepancy 
between X-ray and lensing estimates assuming spherical symmetry, then a simultaneous 
fit of both data sets using a spherical model is not appropriate. 

The non parametric method based on Abel's integral discussed in the previous 
Section has been applied to some clusters in order to stud y the temp e rature and 



density profiles of th e ICM, assum ing spherica l sym metry: Yuan et al.l (l2008l) and 
Kitavama et all l|2004h for RXJ 1347: lNord et all l|2009l ) for Abell 2163; and lBasu et"aLl 
(|2010l ) for A"bell 2204. 
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Recently. iMahdavi and Chand ( 201ll ') derived a model-independent expression for 



the minimum line of sight extent of the hot plasma in a cluster of galaxies. No a priori 
assumptions regarding equilibrium or geometry are required, and the inputs are X- 
ray and SZ data. They applied this method to the Bullet Cluster, and constrained a 
minimum line of sight / plane of the sky axial ratio of ~ 1. 

In deprojecting, assumptions regarding the three-dimensional structure of the DM 
or the ICM are often required, as well as the assumption regarding hydrostatic equilib- 
rium. Although these assumptions could be justified in relaxed and uni-modal clus- 
ters (i.e. clusters for which the mass distribution can be described using a single 
mass clump), they are likely to be too strong for the vast majority of the systems. 
Indeed, numerical simulations often show that clusters deviate systematically from 
hydrostatic equilibrium, especially at large radii {R > R500). The main interpreta- 
tion for this result is that clusters, being relatively young structures, are dynamically 
active and continue to accrete mass from their outskirts. Bulk motions of the gas 
lead t o non-thermal pressure support and thus alter the state of hydrost atic equilib- 
rium ( Piffaretti and ValdarninillioOSi : iLau et al.llioogi ') . In fact, as shown bv lRasia et al.l 



( 20061 ) . X-ray mass estimates of simulated clusters based on the hydrostatic equilibrium 



assumptions are biased low by 10 — 15%. The bias i s however d e pende nt on the gas 
physics implemented in the simulation. For example, iRasia et al.l (|2012l ) find that the 
bias is sensitive to the model assumed to include thermal conduction in hydrodynam- 
ical simulations. Observational results by the CCCP and the LoCuss collaborations 
(e.g. iMahdavi et al. 2008; Zhang ot al. 2010), based on the comparison between X-ray 
and lensing masses of large samples of clusters, seem to support the view emerging 
from the simulations: if we belie ve that the lensing masses are nearly un-biased on 
average (|Meneghetti et a"l]|2010ah . the ratio between X-ray and lensing masses should 
give an indication of the lack of hydrostatic equilibrium in these systems. From the 
above mentioned observational projects, it emerges that X-ray masses are generally 
smaller than the lensing masses by ~ 15%. 



3 Priors 

Besides combining different observational data sets and performing a simultaneous fit 
of a triaxial mass model, it is also possible to consider well motivated priors in order to 
narrow down the parameter space that will be explored when performing the fit. We 
refer the reader to Corless a nd King. (,20081 and Corlcss ct al. (2009) for a thorough 
discussion. Use of priors applies when the fit is performed on a single given data set 
or when the fit is applied to complementary data sets. Of course, the fewer data sets 
one is using and/or the poorer the quality of these data sets, the more one relies on 
priors. Priors sometimes rely on the results from N-body simulations which reasonably 
capture the physics of the DM, being dominated by the gravity and hence relatively 
simple. However, predictions are coming from dissipationless N-body simulations, and 
the physics of baryons may modify predictions, specially in the inner part of galaxy 
clusters where m ost of the baryons are fo und. 

For example. ICorless and Kin3 ( 20081 ) presented a Bayesian MCMC (Monte Carlo 



Markov Chain) method to fit fully triaxial NFW haloes to weak lensing data. Their 
method allows to combine weak lensing data with prior probability functions on the 
model haloe parameters to return parameter and error estimates that reflect the true 
uncertainties of the problem. 
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A similar approach has bee n proposed by Mauro Sereno and collaborators ( Sereno and Umetsul 



I2OI1I : [Sereno and Zitrinllioid ). Using weak and strong lensing data sets, they depro- 
jected the surface density maps to infer the triaxial structure of the cluster. Priors 
were propagated through the inversion by means of the Bayes theorem. The method 
proceeds as follows. As a first step, they obtain the surface mass density of a galaxy 
cluster by strong- lensing modelling of multiple images or/and weak lensing analysis of 
shear and magnification. As a second step, a projected ellipsoidal NFW haloe with ar- 
bitrary orientation is fitted to the convergence map. Finally, the measured distributions 
of projected parameters are recovered using Bayesian statistics. 

Given the measured converg ence map fcobsi the weak lensing XwL function can be 
expressed as jOg uri et al.| |2005 |), 



XWL = t'^^bsCri) - k(ri)] {V-')^ . [fcobs(rj) - k(rj)] (5) 



where V is the inverse of the pixel-pixel covariance matrix. The corresponding like- 
lihood is £wL « exp(-XwL/2)- 

As far as strong lensing is concerned, one can use either a parametric or a non- 
parametric approach in order to retrieve the surface mass density map. 

The final likelihood to be employed in the Bayes' theorem is then £(«;s, rsp , e, Se), 
where Ks,rsp,e,0£ corresponds to the lensing strength, the projected length scale, the 
projected ellipticity and the ellipticity in the plane of the sky of an ellipsoidal NFW 
haloe respectively. Each projected parameter is on turn a function of the intrinsic 
shape and orientation. According to the relevant data-set, the likelihood is £wLi '^SL 
or -Cwl X -CgL for weak, strong or combined strong plus weak lensing analyses, 

respectively. 

We emphasise that strong lensing clusters are typically more elliptical in the core 
than in the weak lensing probed region (see, e. g iMeneehetti et al.ll2010al ). This may 
have important consequences when the likelihoods are combined. 



3.1 Priors on the Axis Ratio 

Using a spherical mass model is equivalent to putting (5-function priors on both axis 
ratios, which is a strong prior. The opposite would be to impose a weak fiat prior on 
both axis ratios. 

A prior on the axis ratio of a triaxial model has been used in recent works by M auro 
Sereno and collaborators. It is detailed below. We note that lOeuri et al.l l|2005l ) were 
the first to use this prior. 

The distribution of minor to ma jor axis ratios (r^oM.a) can be approximated as 
I Jing and Sutoll2002l : iLee et al.ll2005l ). 



P{VDM,a) « exp 



{VBM,a - VtJ./r 



(6) 



where the parameters of these distributions were obtained from numerical simulations: 
r\it, = 0.54, as = 0.113 and 



N0.07r2M(z)° 



(7) 
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with A/* the characteristic nonlinear mass at redshift z and M^-„ the virial mass. 
The conditional probability for intermediate to major axis ratio 7;dm,6 goes as 



p('7DM,a/'^I5A/,6l'7DM,a) = 



2(1 



2?7DM,a/T?DM,6 - 1 - 7"n 



(8) 



for J7DM,a/'7DM,b ^ '"min = "^^ViTiyv.a^ 0.5], whereas is null otherwise. The distribution 
of axial ratios of the lensing population mimics that of the total cluster population 
( Hennawi erallboOTl ). 

As prior for the intrinsic shape, we consider a flat distribution for the axial ra- 
tios in the range ryniin < 'yoM.a ^ 1 a^nd 77DM,a ^ '7dm. fc ^ 1- Probabilities are de- 
fined such that the marginalised probability P(??DM,a) smd the conditional probability 
^('7DM,&i'7DM,a) s^rs constant. The probabilities can then be expressed as 

P('?DM,a) = 1/(1 - '^min) (9) 

for the full range rymin < '7dm, a ^ 1 

P('7DM,6l'7DM,a) = (1 " '?DM,a)~"^ (10) 
for '7dm, 6 ^ 'yDM.a smd zero otherwise. The resulting probability for '7dm, 6 is then 

p('7DM,fc) = In f ) . (11) 

A flat distribution allows also for very triaxial clusters (77DM,a ^ '7DM,f) ^ l)i which 
are preferentially excluded by A'^-body simulations. Therefore, 



3.2 Strong Lensing Clusters 



Different authors, based on simulations (iHenna wi et ablboOTl : ICorless and Kin^l2007l : 
lOeuri and BlandfordlboogUMeneehettT et al .l2010ah . investigated to whi ch degree strong 



lensing (SL) clusters constitute a biased population of galaxy clusters. iHennawi et al 



(fioOT 'l found that strong lensing clusters have 3D concentrations 18% higher than 
the typical cluster with similar mass. Besides, strong lensing clusters are found to be 
triaxial and viewed preferentially along their major axis. Therefore, we expect an ad- 
ditional bias to exist in the dist ribution of 2D concen trations (the quantity to which 
lensing is sensitive to). Indeed, IHennawi et al.l (|2007l ) found that they have 2D con- 
centrations which are 34% higher than the typical cluster. The bias i n concentration 
of strong lensing clusters may be even higher than the 18% found by IHennawi et al.l 
( 20071 ). iMeneghetti et al] 1 2010ah showed that the strongest gravitational lenses (i.e. 



characterised by large lensing cross sections) are typically affected by much larger con- 
centration bias (up to 100%). The bias also depends on the cluster redshift, being 
stronger at those redshifts that are least favourable for strong lensing. Indeed, at these 
redshifts, only clusters very elongated along the line of sight seem able to produce 
strong lensing events. 

It is now well established that strong lensing clusters constitute a biased population 
of triaxial haloes whose major axis is preferentially aligned with the line of sight, 
boosting the lensing efficiency. Therefore, when studying galaxy clusters presenting 
strong lensing features, using a prior on the angle between its major axis and the line 
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of sight is reasonable and has been used in a numb er of studies (e.o. lCorless et aLll2009l : 
ISereno and Umetsullioill : ISereno and Zitrin|[20li ') . 

If 6 represents the angle between the major axis of the haloe and our line of sight, 
this orientation bias can be expressed as: 



p(cos 6) oc exp 



(cose - 1)" 
2ai 



(12) 



A value of ag = 0.115 can be repres entative of the orientation bias for massive strong 



lensing clusters i Corless et al.l [20091 ) 



Apart from this strong orientation bias, it is worth noting that strong lensing clus- 
ters do not exhibit a significant excess of triaxialit y and display nearly the same di stri- 



bution of axis ratios as the total cluster population (|Hennawi et al.-2007. : .Menee:h c tti et al.l 
l2010al ). 

Having mentioned the existence of this orientation bias, it is worth noting that SL 
clusters do not constitute an homogeneous population: all SL clusters are not triaxial 
haloes with their 

major axes pointing right at us. Actually, observing deep enough at any X-ray 
lu minous cluster, it i s very likely that SL features will be detected. As mentionned 
by iMeneghetti et al.l ( 2010ah . the orientation bias is a growing function of the lensing 



cross section. Large Einstein radii clusters will tend to be more biased in their ori- 
entation compared to smaller Ei nstein radii clus t ers. Fo r example, in the sample of 7 
strong lensing clusters studied bv lNewman et al.1 (|2012al 'l. only 1 appears to present an 



Another way to 'boost' the lensing efficiency is through merger. Indeed, mergers 
provide an efficient mechanism to substantially increase the strong lensing efficiency 
of individual clusters (Zitrin et al. 2012; R cdlich et a l. 2012). Therefore, both merging 
and elongated clusters along the line of sight are overrepresented in strong lensing 
cluster samples. In most cases, the optical properties should be able to disentangle 
both populations. 



3.3 Priors from the Mass-Concentration Relation 



A/^-body siinulations dOeuri and BlandfordI l2009l : iMaccio et aP I2OO8I : ICao et all I2OO8I : 
iDuffy et al.|[2008l : |Prada et al.ll2012l ) have provided a picture of the expected properties 
of dark matter haloes. Results may depend on parameters such as the overall normal- 
isation of the power spectrum, the mass resolution, and the simulation volume. The 
dependence of haloe concentration on mass and redshift can be adequately described 
by a power law. 



c = AiM/Mpi^otf{l + zr- (13) 

As reference, we follow [Puffv et al] 1 20081 ). who used the cosmological parameters from 
WMAP5 and found {A, B, C} = {5.71±0.12, -0.084±0.006, -0.47±0.04} for a pivotal 
mass Mpivot = 2 x 10^'^ Mq /h in the redshift range — 2 for their full sample of clusters. 
The scatter in the concentration about the median c(M) relation is lognormal, 



p{lnc\M) 



2-K 



exp 



Inc- lnc(A/) 



(14) 
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with a di spersion gflo R m C2r) o) = 0.15 for a full sample of clusters ( Duffy et al.ll200^ 1. 
Recently, IPrada et al.l (|2012l 'l claimed that the dependence of concentration on haloe 
mass and its evolution can be obtained from the root-mean-square fluctuation ampli- 
tude of the linear density field. They noticed a flattening and upturn of the relation with 
increasing mass and estimated concentrations for galaxy clusters substantial ly larger 
than results reported in Eq. (|13|l . However, more recently. iLudlow et al studied 
how the dynamical state of dark matter haloes affects the relation between mass and 
concentration. When considering only dynamically relaxed haloes, they find that the 
aforementioned upturn disappears. Finall y, Meneghett i fc Ra sia (submitted) show that 
the high amplitude and the upturn of the lPrada et al mass-concentration rela- 

tion can be explained in terms of: i) the different method (compared to other works in 
the literature) to measure the concentration (from the measured circular velocity) and 
ii) the different selection applied to haloes for building the mass-concentration relation 
(clusters are selected by maximum circular velocity). 

In the literature, there is no consensus on the evolution of the c(M) relation with 
redshift, in particular for massive haloes. While several a uthors predict a strong redshift 
evolu ti on of the concentra tions at all mass scales (see, e.o. lBuUock et al ]l200ll : lEke et al.l 
l200ll ). IZhao et al.1 (jiooi) find that the evolution of the concentration of individual 
haloes is not just a function of re dshift but is tightly connected to their mass growth 
rate (see also lWechsler et al.ll2002l l. In particular, the faster the mass grows, the slower 
the concentration increases. Since most of the massive haloes are in a fast mass accretion 
phase at high redshift, they find that the cluster c( M) relation has a very s low re dshift 
evolution. Similar results were found recently by iMufioz-Cuartas et al 

] (|201lh . who 

confirm that the growth rate of the concentration depends on the haloe mass, with 
low-mass haloes experiencing a faster concentration evolution. These authors also find 
that the evolution of the c(M) relation is faster at lower redshifts than at higher 
redshifts. 



3.4 On the Choice of Priors 

The choice of priors used on the parameters of a triaxial model is very important 
because of the inherently under-constrained nature of the problem. This choice should 
be carefully made with respect to the particular problem at hand. For example, if 
one wants to test the mass-concentration relation predicted by the ylCDM scenario, 
adopting a prior based on that relation would not be relevant. On the other hand, 
if one aims to model a sample of dark matter haloes to calculate a mass function, 
a loose mass-concentration prior may be appropriate in order to take into account 
exi sting knowled g e of the cluster and group population into the model. The work 
bv ICorless et al.l ( 20091 ) further illustrates this point: their weak lensing analysis of 



Abell 2204 is able to constrain the ellipticity in the plane of the sky, which is found 
to be slightly larger than what is found statistically in /ICDM simulations. Therefore, 
imposing a prior on the axis ratio derived from simulations for this cluster suppresses 
real information and is not appropriate. This highlights the fact that a prior may be 
well adapted to determine the statistics of a large population of dark matter haloes, 
but it may not be relevant for an individual cluster. 

When considering priors within a Bayesian framework, one can compute the Bayesian 
evidence with or without the prior and see which hypothesis performs better. It is worth 
noting that, the better the data, the lesser the impact of priors in Bayesian methods. 
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4 Three-Dimensional Structure of Galaxy Clusters: Triaxiality 

We present in this Section a general parametric framework intended to simultaneously 
fit X-ray, SZ, and gravitational lensing (both weak and strong) data sets. It is based on 
the works published by An drea Morandi and collaborators. More detail can be found 
in the relevant publications ( Morandi et al.l[2010l . l2011all2011bl : lMorandi and Limousiiil 



l2012l : lMorandi et al.|[2012bl ). 

The lensing effects and the X-ray/SZ emission both depend on the properties 
of the DM gravitational potential well, the former being a direct probe of the two- 
dimensional mass map via the lensing equation and the latter an indirect proxy of 
the three-dimensional mass profile through the hydrostatic equilibrium (HE) equation 
applied to the gas temperature and density. In order to infer the model parameters of 
both the IC gas and of the underlying DM density profile, we perform a joint analysis 
of lensing and X-ray/SZ data. We briefly outline the methodology in order to infer 
physical properties in triaxial galaxy clusters: 

— We start with a general ised Navarro, Frenk a nd White (gNFW) triaxial model of 
the DM as described in ljing and Sutj (|2002l ). which is representative of the total 
underlying mass distribution and depends on a few parameters to be determined, 
namely the concentration parameter C200, the scale radius _Rs, the inner slope of 
the DM 7 , the two axis ratios (r?nM a and 7?DM,b) and the Euler angles ^p, 9 and 4> 

— following iLee and Sutol (|2003l . bood l. we recover the gravitational potential and 
two-dimensional surface mass S (Equation I32p of a dark haloe with such triaxial 
density profile 

— we solve the generalised HE equation, i.e. including the non-thermal pressure Pnt 
(Equation I26|l . for the density of the IC gas sitting in the gravitational poten- 
tial well previously calculated, in order to infer the theoretical three-dimensional 
temperature profile T 

— we calculate the SZ temperature decrement map AT{p) (Equation I30|) and the 
surface brightness map S'x (Equation [29]) related to the triaxial ICM haloe 

— the joint comparisons of T with the observed temperature, of with the observed 
brightness image, of AT{u) with the observed SZ temperature decrement, and of S 
with the observed two-dimensional mass map give us the parameters of the triaxial 
ICM and DM density model. 



4.1 ICM & DM Triaxial Haloes 

We describe the DM and ICM as ellipsoids oriented in an arbitrary direction on the 
sky. We introduce two Cartesian coordinate systems, x = {x, y, z) and x' = (x' , y' , z'), 
which represent respectively the principal coordinate system of the triaxial dark haloe 
and the observer's coordinate system, with the origins set at the centre of the haloe. 
We assume that the z'-axis points along the line of sight direction of the observer and 
that the x' , y axes identify the directions of West and North, respectively, on the plane 
of the sky. We also assume that the a;, y, z-axes point along the minor, intermediate 
and major axes, respectively, of the DM haloe. We define ■(/>, 9 and (/> as the rotation 
angles about the a;, y and z axis, respectively (see Figure[T]). Then the relation between 
the two coordinate systems can be expressed in terms of the rotation matrix M as 



x' = Mx, 



(15) 
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Fig. 1 The orientations of the coordinate systems. The Cartesian axes {x, y, z) represent 
the DM and ICM haloe principal coordinate system while the axes {x',y',z') represent the 
observer's coordinate system with the z'-axis aligned with the line of sight (l.o.s.) direction. 
We define ip, 8 and <f) as the rotation angles about the x, y and z axes, respectively. The labels 
N and W indicate the position of the North and West, respectively, on the plane of the sky. 



where # denotes matrix multiplication and M represents the orthogonal matrix corre- 
sponding to counter-clockwise/right-handed rotations Mx{tp), My (6), Afz(0) with Euler 
angles il},9,(j), and it is given by: 



M = M:,{i,)#My{e)#M^{(t>) , 



(16) 



where 



M^(V') = 



My {9) 



1 

cos i/; — sin ijj 
sini/; cosijj 

cos 9 sin 61 

1 
— sin 9 cos 9 



(17) 



(18) 



M40) = 



cos cf) — sin cj) 
sin (j) cos (j) 
1 



(19) 



Figure [T] represents the relative orientation between the observer's coordinate sys- 
tem and the haloe principal coordinate system. 

In order to parametrise the cluster mass distrib ution, we consider a triaxial gener- 
alised Navarro, Frenk & White model (gNFW, e.g. Ijing and Sutol[20oi ): 
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where -Rs is the scale radius, 5c is the dimensionless characteristic density contrast with 
respect to the critical density of the Universe pc,z at the redshift z of the cluster, and 
7 represents the inner slope of the density profile; pc,z = 3H{z)'^ /8tvG is the critical 

1/2 

density of the universe at redshift z, Hz = Ez Hq, Ez = [i7M(l + + > 

_ 200 e^oo 



where C200 = R20o/Rs is the concentration parameter, with ( Wvithe et alJlioOlh : 



F(y,7)= r s'-^l + sy-^ds. (22) 







The radius R can be regarded as the major axis length of the iso-density surfaces: 
i?2 ^ f 4 + 17 + 4V {a<b<c). (23) 



We have previously defined »7DM,a = Oj/c and ?7DM,f) = b/c as the minor-major and 
intermediate-major axis ratios of the DM haloe, respectively. 

The gravitational potential of a dark haloe with the triaxial density profile (Equa- 



tion 20)1 can be written in terms of complex implicit integrals I Binnev and Tremainel 



119871 ). While num erical integration is r equired in general to obtain the triaxial gravi- 
tational potential, iLee and Sutol (|2003l ) retrieved the following approximation (which 
holds for small eccentricities for the gravitational potential $ under the assumption of 
triaxial gNFW model for the DM (Equation [20j : 

^^(u) Co F^iu) + Co ^^F,iu)+Co i^ll^Eq±A^ p,^^) , (24) 

with u = v/ Rs, Co = 47rG(5ePc(^)fis 1 and t he three functions, Fi(ii), F2(u), and F^iu) 
have been defined in lMorandi et all ( 2010l ). (ef,) and Ec (ec) are the eccentricity of 



DM (IC gas) with respect to the major axis (e.g. ef, = \/l — (fe/c)^). 

The work of lLee and Sutol (|2003l ) showed that the iso-potential surfaces of the triax- 



ial dark haloe are well approximated by a sequence of concentric triaxial distributions 
of radius -Ricm with different eccentricity ratio. For -Ricm a similar definition as for R 
holds (Equation I23|) . but with eccentricities and tc. Note that — ei,{ei,,u,y) and 
Ec = ec(ec,u,7), unlike the constant 6^,60 for the adopted DM haloe profile. In the 
whole range of u, thl^h [^c/ec) is less than unity (~ 0.7 at the centre), i.e ., the intraclus- 
ter ga s is altogether more spherical than the underlying DM haloe (see lMorandi et al] 
(|2O10l ) for further details). 

The iso-potential surfaces of the triaxial dark haloe coincide also with the iso- 
density (pressure, temperature) surfaces of the intracluster gas. Thi s is simply a direct 
consequence of the X-ray shape theorem ( Buote and Canizarej[l993 ): the HE equation 



(gni yields 

VP X V<5 = Vpgas X V<P = 0. (25) 
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4.2 X-ray, SZ and lensing equations 



For th e X-ray analysis we rely on a generalisation of the HE equation ( Morandi et al.l 
which accounts for the non-thermal pressure Pnt and reads: 

VPtot = -/5gasV<? (26) 

where pgas is the gas mass density, <1> is the gravitational potential, Ptot = Ah + Ait- 
We implemented a model where Pnt is a fractio n of the total pres sure Ptot , and we set 
this fraction to be a power law with the radius ( Shaw et al.ll201ol '): 



^=C(i?/i?20o)" . (27) 

Ptot 

Note that X-ray data probe only the thermal component of the gas Pjjj = rie kT, 
k being the Boltzmann constant. From Equations (|26|l and (|27p we point out that 
neglecting Pnt (i.e. Ptot = Pth) systematically biases low the determination of cluster 
mass profiles. This effect increases at larger radii, where the contribution of the gas 
motion is larger. 

Given that Equation (|26p is a first order differential equation, we need a boundary 
condition on the pressure, P, which represents the pressure at P20O1 and it is an 
unknown parameter to be determined. 

To model the electron density profile in the triaxial ICM haloe, we use the following 
fitting function , which corresponds to a simplified version of the function given by 
IVikhlinin eTsI\ (|2006al '): 



ne(Picm) = no (PicmKJ-'(l + Pfcm Ac J (1 + ^fcm/'^l ) ""^^ (28) 

with parameters {uq, Tqi ,£,S, rc^ ,v). We computed the theoretical three-dimensional 
temperature T by numerically integrating the equation of the HE (Equation [26]), as- 
suming triaxial geometry and a functional form of the gas density given by Equation 

The observed X-Ray surface brightness Sx is given by: 

= 4^(l + ,)4 -^(^P*-°j'^) / "^"P ' (29) 

where yl(rp;.oj, Z) is the cooling function. Since the projection on the sky of the plasma 
emissivity gives the X-ray surface brightness, the latter can be geometrically fitted with 
the model ne(Picm) of the assumed distribution of the electron density (Equation [28]l 
by applying Equation l|29p . This has been accomplished via simulated Chandra spectra, 
where the current model is folded through response curves (ARF and RMF) and then 
added to a background file, and with absorption, temperature and metallicity measured 
in that neighbouring ring in the spectral analysis. In order to calculate A{T*j.^j, Z), we 
adopted a MEKAL model for the emissivity. 

The thermal SZ effect is expressed as a small variation in the temperature AT{u) 
of the CMB as a function of the observation frequency: 



AT{u) _ GT 



Pe{v)f{v-T{v))dz' (30) 
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where aj' is the Thomson cross-section, Pe(i') = ne(r)kre(r) is the pressure of the 
electrons of the ICM at the volume element of coordinate r, k is the Boltzmann con- 
stant, and Tcmb ~ 

2.725 K. 

f{i';T{r)) takes into account the spectral shape of the SZ effect and it reads; 

/(z.; r(r)) = {^^^ - 4)(1 + Of{x; T)) , (31) 

where x = hv/kT^^^^ accounts for the frequency depen dence of the SZ effect, and for 
the relativistic corrections related to the term Of{x,T) I Itoh et allll998l ). Note that in 



Equation pO|) we account for the implicit dependence of f(y;T{v)) on the radius. 

Next, the two-dimensional SZ model AT{u) is convolved with the instrumental 
point-spread function and the measured transfer function. In practise, the transfer 
function convolution is performed via multiplication in the Fourier domain. This fil- 
tering significantly reduces the peak decrement of the cluster and creates a ring of 
positive flux at r ~ 2 arcmin. This filtered model is then compared to the observed SZ 
temperature decrement map. We also calculated the noise covariance matrix C among 
all the pixels of the observed SZ temperature decrement map through 1000 jackknife 
realisations of our cluster noise. In this perspective we assumed that the noise covari- 
ance matrix for SZ data is diagonal, as this was shown to be a good assumption in 
ISavers et al] (|2011al ). 



For the lensing analysis the two-dimensional surface mass density S can be ex- 
pressed as: 

/>oo 

p{R)dz' (32) 

We also calculate d the covarianc e matrix C among all the pixels of the reconstructed 
surface mass fsee lMorandi et al] |2011b ) for further details). 



4.3 Joint X-ray-|-SZ-|-Lensing Analysis 

The probability distribution function of model parameters has been evaluated via 
Markov Chain Monte Carlo (MCMC) algorithm, by using as proposal density a likeli- 
hood C and a standard method for rejecting proposed moves. This allows to compare 
observations and predictions, and to infer the desired physical parameters. The like- 
lihood has been constructed by performing a joint analysis for lensing and X-ray/SZ 
data. More specifically, the system of equations we simultaneously rely on in our joint 
X-ray -|- SZ -|- lensing analysis is: 



r(c2oo, -Rs, 7, jyoM^a, '7dm,6. V', ^, 0, ^^0, ''ci , £, 5, fcj, V, P) 
'S'x(c200, -Rs, 7i mM,a, VDM.b, "0, ^o, , £, 5, , u) 
^r(c200, -Rs, 7, VDM,a,VDM,b, V', <t>, "•0, ''ci , £, 5, rc2,V, 5, n, P) 

S(C200, -Rs, 7, VDM,a,VDM,b, V', 0) (33) 

where the parameters C200 (concentration parameter), Rs (scale radius), 7 (inner DM 
slope), »7DM,a (minor-major axis ratio), 77dm, & (intermediate- major axis ratio), and 
■4),9,<j) (Euler angles) refer to the triaxial DM haloe (Equation I20p : the parameters 
no, rci , e, (5, , w refer to the IC gas density (Equation [28]); (normalisation and 
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slope, respectively) refer to the non-thermal pressure (Equation (27]); and P to the pres- 
sure at -R20O1 which is a boundary condition of the generalised HE equation (Equation 
[HI. 

In this triaxial joint analysis the three-dimensional model temperature T is recov- 
ered by solving equation (|26p and constrained by the observed temperature profile; the 
surface brightness is recovered via projection of the gas density model (Equation I29[l 
and constrained by the observed brightness; the SZ signal is deduced via projection 
of the three-dimensional pressure (Equation I30|l and constrained by the observed SZ 
temperature decrement; and the model two-dimensional mass density S is recovered 
via Equation (I32|) and constrained by the observed surface mass density. 

Hence the likelihood £ oc exp(— x^/2), and reads: 



X 



:^ = xix + Xx,S + Xsz + xfcns (34) 

O O Q Q O 

with Xx T' Xx S' Xsz ^'^'i Xions being the x coming from the X-ray temperature, X- 
ray brightness, SZ temperature decrement and lensing data, respectively. We note that, 



when both weak (lensing WL) and strong lensing (SL) data are available (jMorandi et al.l 

|2011b|), xLs = XwL+XsL- 

For the spectral analysis, Xx T equal to: 

X^T = E (35) 

7 — 1 ^ ■ ■ 



Tpioj,i the azimuthally-averaged projection (following IMazzotta et ani2004l ') of the the- 
oretical three-dimensional temperature T; the latter is the result of solving the HE 



equation, with the gas density ne(-Ricm). 
For the X-ray brightness, Xx S reads: 



7 i=l -^i* 



with 5x,i and S"^ j theoretical and observed counts in the ith pixel of the jth image. 
Given that the number of counts in each bin might be small (< 5), then we cannot 
assume that the Poisson distribution from which the counts are sampled has a nearly 
Gaussian shape. The standard devi ation (i.e . , the square-root of the variance) for this 
low-count case has been derived bv lGehreld l|l986h : 



^S,^ = 1 + ^ S* + 0.75 (37) 



which has been proved to be accurate to approximately one percent. Note that we 
added background to Sx,i as measured locally in the brightness images, and that the 
vignetting has been removed in the observed brightness images. 
For the SZ (lensing) constraint D, the Xd contribution is: 

Xd = [D-D*]'C"^[D-D*] , (38) 



where C is the covariance matrix of the two-dimensional SZ temperature decrement 
(projected mass density), D* are the observed measurements of the two-dimensional 



19 



SZ temperature decrement (projected mass density) in the ith pixel, and D is the 
theoretical 2D model. 

Errors on the individual parameters have been evaluated by considering average 
value and standard deviation on the marginal probability distributions of the same 
parameters. 

So we can determine the physical parameters of the cluster, for example the 3D 
temperature T, the shape of DM and ICM, just by relying on the generalised HE equa- 
tion and on the robust results of the hydrodynamical simulations of the DM profiles. 
In Fig. [5]we present an example of a joint analysis for T, Sx, ^T{i/) and S: for Sx, 
AT{i') and S the ID profile has been presented only for visualisation purpose, the 
fit being applied on the 2D X-ray brightness/SZ/surface mass data. Note that in the 
joint analysis both X-ray, SZ and lensing data are well fitted by our model, with a 
X?cd = 1-04. 



4.4 Limitations 



The BCG: The centre of galaxy clusters is usually populated by a bright central galaxy 
(BCG). Actually, the BCG has an influence on the formation of multiples images 
( Meneghetti et a"l]|2003l : iDonnarumma et al.ll201ll ). and the physical proc esses taking 



place in the BCG substantially influence the X-ray gas ( Gitti et al.ll2012r i. In the tri- 
axial framework described in this Section, we have removed the central 25 kpc of the 
data in the joint analysis, to avoid the contamination from the BCG. However, mod- 
elling properly the BCG contribution is essential in order to probe the dark matter 
distribution in the very centre. 

In galaxy cluster Abell 1703, the strong lensing analysis bv lLimousin et ahl (|2008h 
takes into account the stellar contribution of the BCG in a parametric mass modelling 
aimed at constraining the underlying smooth dark matter component distribution, in 
particular the inner slope of the dark matter distribution. 

A more advanced approach which requires high quality spectro scopic data is to 
comb i ne lensing observations with the stellar ki nematics of the BCG (|Sand et al.ll200 j . 
|2004 I2OO8I : iNewman et al.1 lioill . I2012all2012bh . These studies found density proflles 
shallower than canonical NEW models at radii < 30 kpc, comparable to the effective 
radii of the BCG. 



Numerical simulations bv iDubinskil ( 19981 ) suggested that the measured velocity 



dispersion proflle strongly depends on the line of sight, where the central value peaks 
between 300 and 450kms~^ depending on the considered line of sight. However, recent 
observatio ns do not confirm this cl aim. In the sample of 7 massive strong lensing clusters 
studied bv lNewman et all (|2012al ). the observed velocity profiles display a very homo- 
geneous shape which are mutually consistent . Similar conclusions is reached within 
the SAURON project (|Cappellari et al.ll2007h . which found that in giant slowly /non 



rotating ellipticals (including BCGs), the average shape of the velocity ellipsoid does 
not differ by more than ~ 10 per cent from a spherical shape. 

The gas mass component: When studying the dark matter distribution, we did not 
subtract the mass contribution from the X-ray gas to the total mass. Nevertheless, the 
contribution of the gas to the total matter is small: the measured gas fraction is 0.06- 
0.07 in the spatial range 30-400 kpc, and the slope of the density profile is very similar 
to that of the DM beyond a characteristic scale ~ 20-30 kpc, a self-similar property of 
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Fig. 2 Example of the joint analysis for T, Sx, ^T{u) and S, from iMorandi et alj l l2012bl '). 
In the upper panel we display the two quantities which enter in the X-ray analysis (Equation 
I35I I: the observed spectral projected temperature T*^^. ^ (big points with errorbars) and the 
theoretical projected temperature Tpi.oj,m (diamonds). We also show the theoretical 3D tem- 
perature (solid line), which generates Tp^oj ^ through convenient projection techniques. In the 
second panel from the top we display the two quantities which enter in the X-ray brightness 
analysis (Equation I36I I: the observed surface brightness profile 5^ (points with errorbars) and 
the theoretical one 5x (solid line). In the third panel from the top we display the two quan- 
tities which enter in the SZ temperature decrement analysis (Equation [SOj : the observed SZ 
temperature decrement profile (points with errorbars) and the theoretical one AT{i') (solid 
line). Both observed and theoretical SZ temperature decrement are convolved with the transfer 
function: note that this filtering significantly reduces the peak decrement of the cluster and 
creates a ring of positive flux at r ~ 2 arcmin. In the lowest panel we display the two quanti- 
ties which enter in the lensing analysis (Equation JSSjl): the observed surface mass profile E* 
(points with error bars) and the theoretical one S (solid line). Note that for surface brightness 
(surface mass) and the SZ data the ID profile is presented for visualisation purpose only, the 
fit being applied on the 2D data. Moreover, for the surface brightness we plotted data referring 
to the observation ID 6880. The virial radius corresponds to a scale length on the plane of the 
sky of ~ ?7DM,a ■ H200 ~ 2240 kpc. 
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Fig. 3 Degeneracies between 9 and ^. 6 represents the angle between the major axis of the 
haloe and our line of sight, and ^ is related to the amount of non-thermal component consid- 
ered in the model. The solid (dashed) line represent the 1 (2)-(t error region, while the point 
represents the best fit value. 



the gas common to cool core clusters i Morandi and Ettorill20o"7l l. suggesting that the 
assumption to model the tota l mass as a gNFW is reliable. Similar conclusion s have 
been reached bv lBradac et al.l ( 20081 ) and ISommer-Larsen and Limousin! I 2010| ). 



4.5 Degeneracies & Priors 



Degeneracies arise between the different parameters involved in the modelling. Fig. |3] 
presents the degeneracies expected between 6 and ^. 6 represents the angle between 
the major axis of the haloe and our line of sight, and ^ is related to the amount of non- 
thermal component considered in the model. In Fig. |4]we present the joint probability 
distribution among different parameters in our triaxial model for Abell 1689. 

Regarding the use of priors, we point out that, when combining complementary 
data sets as described in this Section, we do not need to rely on any priors like those 
discussed in Section [S] This is welcome since priors may be potentially biased due to 
our incomplete understanding of cluster physics. 
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5 Three-Dimensional Mass Distribution in A1689 (and Other Clusters) 

We aim to illustrate how spherical and triaxial modelling can lead to different haloe 
parameters. For this purpose, we focus on galaxy cluster Abell 1689 which has been 
extensively studied at different wavelengths. Abell 1689 is a massive galaxy cluster at 
redshift 0.18 with a very large Einstein radius, around 45" for a source redshift ~ 2. 

Before going into detail, we would like to emphasise the complexity of this structure 
an d try to justify why we treat such a complicated structure using a single mass clump. 

iLokas et al.l (|2006h used spectroscopic redshifts to study the kinematics of about 
200 galaxies in the cluster. They showed that the clu ster is probably surrounded by a 
few structures aligned along the line of sight. [Czoskd (2004) reported redshifts for 525 
galaxies, spanning from the centre outward to 3h~^ Mpc. They found only one appar- 
ently distinct group of galaxies that lies about 350 kpc to the northeast of the cluster 
centre. It corresponds to a group of bright galaxies well identified in optical images of 
this cluster. All strong lensing studies have taken into account the gravitational pertur- 
bation it generates, but its contribution to the total mass budget is found to be small. 
The redshift distribution of these galaxies is skewed toward slightly higher redshifts. On 
larger scales (R> lh~^ Mpc), no evidence for any substructures is found: the outskirts 
of Abell 168 9 look rather homogeneo us. However, deep X-ray observations obtained 



with Suzaku I Kawaharada et al.l 2010l l reveal anisotropic gas temperature and entropy 



distributions in the cluster outskirts. 

If Abell 1689 appears as a complex structure, the main mass clump seems to be 
dominant in the mass budget: all strong lensing studies find the mass centre to coincide 
with the brighte st cluster galaxy, w hich also coincides with the peak of the circular 



X-ray emission ( Lemze et al. 20081'). Besides, deep Chandra data have revealed the 



presence of a cool core ( Riemer-Sarensen et al.ir2009l 'l. 



5.1 Spherical Analysis: the Abell 1689 Puzzle 

While it has been proposed as a standard example of a relaxed object in hydrostatic 
equilibrium, Abell 1689 has been the subject of some controversies. To summarise, as 
long as spherical symmetry has been assumed, this cluster has been problematic in two 
ways: 

— High concentration parameters (up to ~ 30) have been derived from lensing analy- 
ses. These large values r epresent a major i nconsistency with the theoretical ylCDM 
expectations (c ~ 3 - 4, iNeto et al.|[2007l : lOuffv et allboosh . 

— The 2D mass derived from X-ray data is only half of the mass found using SL 
estimates. 

In the following, we discuss these two points and then we turn to triaxial models 
and how they have allowed to resolve these issues. 

The concentration parameter of Abell 1689 has been constrained extensively by 
different authors, findin g very different r esults. In weak lensing analyses , first stud- 



ies reported C2oo= 4.8 |King et all 120021 '): 3.5l°-^ jBardeau et al.l l2005h . but these 



analyses likely suffered from dilution of the inner shear profile by foreground clus- 
ter members. On the other hand, very high concentra tions were also infe rred from 



weak lensing analyses b ased on Subaru data: C2oo=30.4 ( Halkola et aklliood ): 22.1^^'^ 



naiysesDc 

I Medezinski et al.ll2007f) . Jsing the same Subaru data but a different algorithm, other 
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auth o rs rep orted a smaller concentration parameter, between 10 and 15 ( Umetsu et alJ 
l2009l . [2011). Weak and stro ng lensin g analyses converged to: C2oo=7.9 ( Clowell2003l 'l: 
7.6+° g jHalkola et al.ir2006l '): 7.6±1.6 ^Limousin et al.lboOTl 'l. values which are still high 
compared to the theoretical expectations. Actually, only a fe w haloes formed in the Mil- 
lennium simulation could reach this value ( Neto et al] |200j ) . These values where foun d 
in agreement with X-ray analyses: C200 = 7.7+2'g ( Andersson and Madeiski|[200i l: 
6.6 ± 0.4 jPeng et aLlliooj ). 

Besides the large variance in the concentration parameter, these studies agree on 
the fact that the mass derived from X-ray meas urement is half of that found from 
str ong gravitational lensing at most radii (see also lLemze et al]|2008l ). 



•ong g 
[Rie i 



lRie mer-S0rcnscn et al. (2009) showed that this discrepancy is reduced if we ex- 
clude a cool clump plus some substructure in the North-Eastern part of the cluster; 
nevertheless a discrepancy still remains in the strong lensing region. 

A way to reconcile the mass derived from X-ray and lensing measurement within 
a spherical mass distri bution is to add the contribution from non-thermal pressure. In 
the case of Abell 1689. lMolnar et al.l (|20ld ') found a contribut ion of about 40%. This is 
larger than the theor etical expecta t ions b y ( Shaw et al.ll201ol ) but consistent with the 
set of simulations bv lMolnar et all ||2010| ). 



5.2 Triaxial Models: Solving the Puzzle 



iQguri et al.l (B ^ , using a triaxial mass model, found that weak lensing measurements 
in Abell 1689 based on Subaru data are indeed compatible with ylCDM if Abell 1689 
repres ents a rare population (~ 6% by number) of cluster-scale haloes. ICorless et al.l 
(I2OO9I ) constrained the triaxial shape of the total mass distribution of Abell 1689 via 
weak lensing data and under a range of Bayesian priors derived from theory, though 
lar ge errors accompany th eir triaxial parameter estimates. 

Morandi et al.l ( 2011ah presented the determination of the intrinsic shapes and the 



physical parameters of both DM and ICM in Abell 1689 by combining X-ray and strong 
lensing data. They showed that Abell 1689 can be described as elongated along the line 
of sight, with a minor-major principal axis ratio equal to 0.42 ± 0.02. They assumed 
that the triaxial ellipsoid is oriented along the line of sight, an assumption justified in 
light of the " orientation bias" of stro ng lensing clusters (see Section r3.2p . A subsequent 
re-analysis of lMorandi et al" I (|2011bl ). by jointly analysing also weak lensing data and 



accounting for the non-thermal pressure of the IC gas, strengthened the view of a 
triaxial cluster elongated along the line of sight, though a bit larger value of the minor- 
major principal axis ratio has been inferred (0.50 ±0.01). Besides, it was shown in this 
work that the large Einstein radius observed was reproduced by a triaxial model with 
ylCDM friendly parameters. In the present review we extend these previous works, by 
allowing the DM and ICM ellipsoids to be oriented in an arbitrary direction on the 
sky (see Section |3] for further details). Our work indicates that Abell 1689 is a triaxial 
galaxy cluster with DM haloe axial ratios r;DM,a = 0.56±0.07 and rjDM b ~ 0.75±0.08, 
C200 ~ 5.27 ± 0.46, and with the major axis slightly inclined with respect to the line of 
sight of 8 — 27.3 ±7.1 deg, in agreement with the predictions of lOguri and BlandfordI 
(|2009h . 



Sereno and Umetsul ( 201 ih developed a method for a three-dimensional analysis of 



the DM haloe via SL and WL data. They re-analysed the wea k lensing convergence 
map of A1689 obtained by lUmetsu and BroadhurstI (120081 ) and lUmetsu et all (|2009l ') 
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Fig. 5 Results of the combined strong and weak lensing analysis. Contour plot of the 
marginalised PDF for M200 and C200 for the dark matter haloe as derived under the prior 
assumptions of either fiat axis ratio distribution and random orientation angles (left panel) 
or Af-body like axis ratio distribution and biased orientation angles (right panel). Contours 
are plotted at fraction values exp(— 2.3/2), exp(— 6.17/2), and exp(— 11.8/2) of the maximum, 
which denote confidence limit regions of 1, 2 and 3(t in a maximum likelihood investigation, 
respectively. The full, long-dashed and dashed lines e nclose the 1, 2 and Sc regions for the pre- 
dicted conditional probabi lity c(M) bv lDuffv et al.l | |2008|) . respectively. Figure adapted from 
ISereno and Umetsul 1 I2OIID . 



on wide-field (~ 30' x 24') Subaru data of A1689. For strong lensing, they employed 
a para metric le nsing analysis method based on the gravlens kernel evaluation of XsL 
values (|Keetonl l2001a,200 Ibl) . The SL likelihood values were computed by comparing 
observed to predicted image positions. Several priors were considered. For the axial 
ratios 7?DM,a and ?7DM,f)i they considered either the A'^-body predictions or a flat dis- 
tribution. For the alignment angle 0, they considered either the biased distribution for 
p{9) (see Eau. ll2|l or a random distribution. For the azimuthal angle they always used a 
random flat distribution. For the mass, they always used a flat prior p(M2oo) = const., 
whereas the a priori PDF for the concentration was fl at in the range < C200 ^ 30 
and null otherwise. The theoretical c{A4) relation from iDufFv et al] l|2008l ) was either 
enforced or neglected. 

Whatever the assumptions on either orientation or shape, the concentration is a bit 
larger but still compatible with theoretical predictions (Fig. [5}. A1689 appears to be a 
quite typical massive cluster with a concentration in agreement with the tail at large 
values of the expected population of clusters of that given mass. Independently of the 
priors, the inferred C200 are only > la away from the predicted median value. Priors 
from A'^-body simulations also help to put an upper bound on the concentrations. 

Mildly triaxial haloes do a better job in fitting data than nearly spherical lenses. 
Values of 0.4 < ?7DM.a ^ 0.8 are more likely than either extremely triaxial or nearly 
spherical shapes. Triaxial shapes predicted by A''-body simulations are in good agree- 
ment with these results. Axial ratios derived assuming a flat distribution are compatible 
at la confidence level with predictions from A''-body simulations. Furthermore, they 
can exclude nearly spherical shapes (t^dm.o ~ '?DM ti ~ 1) a^t the 3a confidence level. 
Finally, they find indications for an orientation bias. 
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Fig. 6 2D mass enclosed within a circular aperture of radius R from lensing data (points with 
error bars), from an X-ray only analysis under the assumption of spherical geometry (solid 
line with the la error grey shaded region), and from a joint X-ray+lensing analysis taking into 
account the 3D geometry (dot-dashed line with the la error cray sh aded region). In this la tter 
case, we see that both estimates agree with each other. (Result from lMorandi et al.ll20l"lal . but 
updated within the full triaxial framework presented in Section |4}. 



The minor-major axis ratio is found in agreement with the findings bv lMorandi et alJ 
confirming a triaxial shape; yet, the value of the concentration parameter 
(c200 = 7.3 ± 0.8) is a bit larger than the values presented in this review (Table [TJ. 

Reg arding the mass discrepancy, the triaxial model proposed by iMorandi et al.l 
(|2011bl) . as well as its extension presented in this review, is able to solve the X- 
ray/lensing mass discrepancy. We plot in Fig. [S]the comparison between the 2D masses 
inferred from X-ray and lensing under spherical and triaxial models. 

We have illustrated that a triaxial mass model for Abell 1689 is able to reconcile 
mass estimates from different probes, as well as to reproduce the large Einstein radius 
using ylCDM friendly parameters. It is worth noting that the value of the inner slope of 
the dark matter density profile, 7, also depends on the adopted geometry. Considering 
Abell 1689, a standard spherical modelling leads to 7 = 1.16 ± 0.04, whereas we find 
7 = 0.92 ± 0.07 using a triaxial mass model (Table [TJ. 



5.3 The ICM Properties in a Triaxial Framework 

The ICM of Abell 1689 within a triaxial model has also been the subject of interest. 
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lOe Filippis et alj (|2005h and ISereno et aD (|2006ll investigated the ICM shape by 
combining X-ray and SZ observations. Given the data sets used (in particular, the 
absence of lensing data), error bars were pretty large and a wide range of geometries 
was possible: they concluded that these data sets are compatible both with a prolate 
and an oblate sha pe. 

More recently. ISereno et al. 1 l|2012l ) implemented some significant improvements by 
using Bayesian inference to determine the intrinsic form: i) The method was still para- 
metric but did not rely anymore on the simple isothermal /3 model for the X-ray data. 
The employed profiles can mimic complex features in either the electronic density or 
the temperature profile for the X-ray data, ii) Instead of the central Compton param- 
eter I/O I they considered the more reliable integrated Comp ton paramet e r, in) Even if 
astronomical deprojection is an under-constrained problem I Serenol[2007l ). ISereno et aP 
( 20121 ) could infer the 3D structure of the cluster using a Bayesian method wit hout 
assuming any speci fic configuration. On the other hand in lOe Filippis et"aL I l|2005h and 



ISereno et al.l (|2006l ). the 3D distribution was assumed to be either triaxial and aligned 
with the line of sight or prolate or oblate. 

This method was applied to Abell 1689, where SZ and X-ray observations cover 
in detail a region R < 1 Mpc. The 3D electron density and temperature were mod- 
elled with parametric profiles. Distributions were assumed to be coaligned and el- 
lipsoidal, with constant ecce ntricity and orienta tion. Intrinsic profiles were taken from 
IVikhlinin et al ] l|2006bl l. and lEttori et~aH (|2009l ). The metallicity was fixed to the mean 
observed value. The elongation enters when 3D profiles are projected into the plane 
of the sky, so that fitting at once X-ray surface brightness, temperature, and the inte- 
grated Compton parameter, one can infer as well as the parameters describing the 
distribution. 

The combined X-ray plus SZ analysis allows to infer the width of the cluster in the 
plane of the sky (parametrised in terms of the ellipticity e) and its size along the line 
of sight (expressed as the elongation e/\). These two observational constraints have to 
be used to infer t he int rinsic shape of the cluster and its orientation. 

ISereno et al.l (|2012h found a minor to major axial ratio for the ICM of 0.7±0.15, 
preferentially elongated along the line of sight. The hydrostatic equilibrium is not in- 
volved in their analysis since they do not derive the mass of t he cluster. This value 
of the axis ratio is in good agreement with the one inferred bv iMorandi et aD ( 2011bl ') 
(assuming generalised hydrostatic equilibrium), who quoted a ratio for the ICM be- 
tween 0.66 and 0.77 (since the method includes a radial variation of the axial ratios). 
This suggests that both approaches are consistent, and that the ICM in Ab ell 1689 is 
not far from being in hydrostatic e quilibrium. Note that Sereno et al ] (|2012l ) use XMM 
data for the temperature, whereas iMorandi et al.l ( 2011a ) use a Chandra temperature, 
which is lower than the XMM temperature by 10-20% due to calibration reasons. 

It is worth noting that several processes like radiative cooling and turbulences are 
affecting the shape of the gas distribution more than the dark matter shape. 



5.4 Overview of the Present Results on Clusters Triaxiality 



So far, four clusters have been studied within the full triaxial framework described in 
SectionH M ACS 1423, Abell 1689 , AbeU 383, and A beU 1835. The resuhs publish ed for 
MACS 1423 (jMorandi et aLlboid ) and AbeU 1689 (jMorandi et al.ll2011all2011bh were 
derived using a triaxial model where the haloe's major axis was aligned with our line of 



28 



Table 1 Best-fit model parameters for tfie four clusters for which a full triaxial modelling 
exists. Error bars correspond to Icr confidence level. The lines 1 — 8 refer to the best fit 
parameters of the DM haloe: C200 (concentration parameter), _Ra (scale radius), 7 (inner DM 
slope), r;DM,a (minor-major axis ratio), ?Jdm 6 (intermediate-major axis ratio), and 1/1,6,0 
(Euler angles). The lines 9 — 14 refer to the best fit parameters no, , £, 5, , of the IC 
gas density, while the lines 15 — 16 to the best fit parameters J, n (normalisation and slope, 
respectively) of the non-thermal pressure. Finally, the last line refers to the best fit parameter 
P of the pressure at -R200 1 which is a boundary condition of the generalised HE equation (see 
relevant Equations in Section |4}. Note that only in the case of Abell 1835 all these parameters 
are constrained. This is due to the inclusion of SZ data. 



Cluster 


Abell 1835 


Abell 383 


Abell 1689 


MACS 1423 


C200 


A QO J- r\ A A 

4.v5Z iL U.44 
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Q IT/ 1 -I r. 
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891.0 ± 114.3 
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7 


1.01 ±0.06 


1.02 ±0.06 


0.92 ± 0.07 
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VDM,a 


0.59 ±0.05 


0.55 ±0.06 


0.56 ± 0.07 


0.62 ±0.04 


'7dm, 6 


0.71 ±0.08 


0.71 ±0.10 


0.75 ± 0.08 


0.72 ±0.06 


1/) (deg) 


3.8 ±4.6 


-13.6 ±5.5 


-35.5 ± 13.7 


-34.4 ±5.4 


e (deg) 


18.3 ±5.2 


21.1 ± 10.1 


27.3 ± 7.1 


34.7 ±8.7 


<f> (deg) 


-55.0 ±6.9 


-16.9 ± 15.9 


-11.1 ±6.7 


-72.3 ±8.3 


no (cm~^) 


0.018 ±0.002 


0.063 ± 0.003 


0.017 ±0.001 


0.15 ±0.02 


rci (kpc) 


117.7 ± 10.1 


26.4 ± 1.7 


119.3 ±5.3 


20.6 ±3.1 


e 


0.68 ±0.02 


0.55 ±0.01 


0.72 ± 0.02 


0.55 ±0.02 


S 


0.82 ±0.03 


0.02 ±0.01 


0.33 ± 0.01 


0.02 ±0.01 


rc2 (kpc) 


1674.3 ± 266.7 








V 


0.44 ± 0.04 










0.177 ±0.065 


0.11 ±0.05 


0.24 ±0.05 


0.08 ± 0.03 


n 


0.77 ±0.21 











P (erg/cm"^) 


(2.7 ±0.7) X 10-13 









sight. Since then, the algorithm was improved and we report the results for MACS 1423 
and Abell 1689 derived using the full triaxial framework. Results of the parameters for 
each cluster are given in Table [T] As discussed in the relevant papers, the parameters 
obtained through a triaxial model can be very different from the one obtained through a 
spherical model. In particular, regarding the concentration parameter, similar to what 
has been found for Abell 1689, lower values are inferred within a triaxial framework. 
Shown on Fig. [7| are some key properties of these four clusters: the inner slope of the 
dark matter density profile, and the masses and concentrations. The comparison with 
results from numerical simulations will be discussed in Section [6. II 



6 Discussion 

6.1 (Limited) Comparison with Simulations 

The parameters of a triaxial model obtained for galaxy clusters are to be compared to 
results from numerical simulations in order to test cosmological models. Such compar- 
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Fig. 7 Properties of the four clusters studied within the full triaxial framework described in 
Section^] Left: inner slope of the dark m atter density profile. Right: Masses and concentrations . 
Results from ylCDM based simulations jNeto et al.ll2007l : lDuffv et al.|[200g : iPrada et al.ll2012fl 
are shown. Error bars correspond to Icr confidence levels. We caution that this comparison can 
be misleading since results from numerical simulations are based on a spherical assumption. 



isons are routinely performed in the spherical case. In order to build a density profile, 
simulators usually do a spherical fit to their elongated profiles: they determine the en- 
closed mass or density assuming spherical symmetry even though the haloes are found 
to be triaxial and asymmetric (Section 1 Since triaxial analyses will have a growing 
importance on the observational side in the future, we advocate the need for simu- 
lations to be analysed within a triaxial framework, allowing reliable and meaningful 
comparisons. 

Actually, the bias induced by fitting a spherical analytic formula to an elongated 
mass profile should be compared to the intrinsic scatter m the concentration of NFW 
haloes. 

Having stressed the possible limitations of comparing the triaxial observational 
results with spherical numerical results, we can go back to Fig. [7| where some key 
properties of four clusters studied within a full triaxial framework are presented: the 
inner slope of the dark matter dens ity profile, and t he masses and concentrations. 
(Sph erical) theoretical expectations ( Neto et al] l2007l : lOuffv et al.l l20oi : IPrada et al.l 
are shown on top of the triaxial results. 

Regarding inner slopes, one should keep in mind that, from the theoretical side, 
numerical simulations did not converge on a proper treatment of the baryonic compo- 
nent which is a difficult task. In this respect, when the baryonic component is taken 
into account, different simulations infer either a steepening or a fiattening of the dark 
matter inner density profile (see the revue on cluster cores in this volume). 

We have mentioned that simulations do not provide the relevant material to com- 
pare triaxial results to theo retical expectation s. Moreover, as far as massive haloes are 
concerned, statistic is poor. lOeuri et al.l (|2009l ) argued that the probability distribution 
of the concentration parameter for very massive haloes (larger that 10^^ M0) and its 
redshift evolution had not been studied with great statistics by N-body simulations. 
It is worth noting that this statement also applies to extreme axis ratios in different 



30 



mass and redshift bins. Therefore, we rely on extrapolations from smaller haloes. More 
reliab le theoretical predictions at the high mass end are needed. Recently. IPrada et al.l 
( 20121 ) argued that previous simulation works underestimated the mean concentration 
at the high mass end and find an upturn in the (c,M) relation. iBhattacharva et al.l 
( 201 ih . in their analysis of an even larger simulation, find no evidence for such an 
upturn but they do find a (c,M) relation that differs in normalisation and shape from 
previous studies that have limited statistics in the upper mass range. Moreover, as 
already mentioned earlier, this upturn may b e due to the inclusion of dynamically 
disturbed merging clusters ( Ludlow et al.l[2012l '). 

The situation is changing and the most recent simulations contain a few hundreds 
of such haloes. For example, lAngulo et al ] [2012 1. based on the Millenium-XXL dark 
matter only simulation, reported 464 haloe s more massive t han 2 x10^^ M©, m agree- 
ment with analytical calculations (see, e. o. iMortonson et al.ll201ll ). 

Another difficulty is that the expected abundance of massive haloes strongly de- 
pends on the algorithm used to find them. lAneulo et all (l2012h show that this abun- 
dance changes by a factor of ~ 2 when Friends-of- Friends (jPavis et al.|[l985l ) selected 
subhaloes and self-bound subhaloes are compared. They argue this is in part a conse- 
quence of large haloes not forming an homogeneous population. Indeed, clusters dis- 
playing similar virial masses present a considerable diversity in shape, concentration, 
and amount of substructure. 

To draw reliable conclusions regarding the properties of clusters and their evolution 
with cosmic time, we would need to study within a triaxial framework a large number 
of (massive) galaxy clusters. These results should be compared to predictions from 
simulations inferred within a triaxial framework. 



6.2 On Which Scales do we Measure Triaxiality? 

The scale on which a triaxial model is constrained depends on the data sets used. 
In the case of Abell 1835, the triaxial model has been constrained up to the virial 
radius thanks to the inclusion of SZ data. This question is linked with the following 
one: out to which radius do we have observational data with a sufficiently high signal 
to noise ratio? Strong lensing is limited to the very central part of a galaxy cluster. 
Regarding X-ray data obtained with current facilities, one can expect to reach Rsoo- 
Regarding weak lensing, beyond R500, the signal to noise ratio becomes small, and 
the lensing signal estimated for a single cluster is likely to be dominated by mass not 
associated with this cluster, either correlated or uncorrelated (Hockstra ct al. 201ll). 
This question is also relevant when we compare with triaxial parameters measured on 
simulated clusters. Ideally, we aim to compare measurements performed on the same 
scales. 



6.3 Implementing Triaxiality in the Mass Distribution or in the Gravitational 
Potential? 

There are two possible approaches when introducing a triaxial model for a galaxy 
cluster: one can either implement triaxiality directly in the mass distribution, or in the 
gravitational potential. 
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An ellipsoidal mass distribution generates isopotential surfaces which are approxi- 
mated by concentric ellipsoids of decreasing axial ratio towards the outer volumes, and 
not of constant shape. Vice versa, if the potential is described by a triaxial ellipsoid, 
then the mass distribution giving birth to this potential will not be described by a 
triaxial ellipsoid. 

Both approaches have their advantages and inconveniences and have been followed 
by different authors. 



6.3.1 A Triaxial Mass Distribution 

We have discussed in the introduction that the DM distribution behaves as an ellipsoid, 
being coUisionless and dominating in term of mass. The dynamics of the gas (which 
is the second largest contribution to the mass after the DM) should be driven by the 
gravitational potential well of the DM. Therefore, it seems well motivated to implement 
triaxiality in the mass distribution as proposed in Section |4l being aware that it has 
its own limitations. 

Recent investigations of gala xy scale dark matter haloes forming in the Aquarius 
simulation ( Springel et al.ll2008l ) show that the shape of a haloe depends strongly on 



its environmen t, the time at wh ich we consider it, and the radius at which we measure 
its shape (Vera; Giro et al.ll2oTlh . At redshift 0, haloes exhibit a variation from prolate 
in the inner regions to triaxial/oblate in the outskirts, which clearly complicates the 
modelling. If these results are derived for galaxy scale dark matter haloes, we can 
expect that some of these complicat ions also arise for cluster scale dark matter haloes. 
Indeed. [Mufioz-Cuartas et al.l ( 201lh found that the asphericity is more pronounced in 



the haloe's central region (i.e. for radii smaller than 30% of the virial radius). However, 
internal regions are more likely to be affected by baryonic physics and therefore the 
inclusion of baryons is needed in order to draw more quantitative conclusions. 

It is worth noting that models which involve a gravitational potential generated by 
an ellipsoidal mass distribution are computationally expensive. 



6.3.2 A Triaxial Potential 



iBuote and Humphrevl (|201lj 2011b^ investigated a different type of model where the 



potential, rather than the underlying mass distribution, is an ellipsoid of constant shape 
and orientation. They inferred analytical formulae for galaxy clusters with triaxial 
potential. These triaxial models themselves lead to a straightforward generalisation of 
analytic spherical models, and hence they are computationally fast. 

In a similar fashion. ISereno et al modelled the gas density distribution as el- 

lipsoid of constant shape and orientation. Note that they did not make any assumption 
about the equilibrium of the gas. 

There is both observational and theoretical work suggesting that this modelling 
has its limitations and which shows that the gas distribution follows isopotential sur- 
faces app roximated by concentric elli psoids of decreasing axial ratio towards the outer 
volumes JSuote and CanizaresI Il994 [l996l: iLee and Sutol booj iHavashi et al.l I2OO7I : 



Morandi et al.ll2O10l : lKawaharall20ia : iDeb et al.ll2012l ). 
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6.3.3 So ? 

It appears that, when the accuracy of the simulations or the observations increases, 
none of the ma ss distribut i on no r the potential can be rigorously described by a triaxial 
model. Indeed, IPeb et all l|2012l ) showed in Abell 1689 that the axial ratio of both the 



DM and the ICM is not strictly constant with radius. 

It is worth noting that, in Abell 1689, Morandi et al. and Sereno et al., using 
different modelling approaches find results for the shape of the ICM in pretty good 
agreement, which proves that both approaches are consistent and deserve additional 
investigation (see comparison in Section (5)1. 



6.4 Large Einstein Radii & the Overconcentration Problem 

The first claim regarding the fact that large Einstein radi i may be a problem for 
the jdCDM scenario was made by iBroadhurst and Barkanal ( 2008al ) . This claim was 



very constructive since it led to a number of interesting studies investigating if strong 
lensing clusters were over-concentrated with respect to theoretical expectations, and 
if the large Einstein radii observed for some extreme clusters were challenging the 
ylCDM scenario. The problems of large Einstein radii and of overconcentrated clusters 
are somehow linked to each other: a very elongated haloe with its major axis aligned 
along the line of sight can lead to a highly concentrated projected surface mass density 
profile resulting in a large tangential critical curve, hence a large observed Einstein 

rad ius. 

lOeuri et all (|2009l ) studied the mass profiles of four clusters by combining strong 



and weak lensing data. The mass profiles were found to be well described by an NFW 
profile. They found values for the concentration parameter that are slightly higher than 
the ylCDM predictions, even after taking into account the lensing bias which includes 
the projection effect. Taking into account the error bars they derived, the value of the 
concentration for each cluster is marginally consistent with theory, and the excess is 
not so strong compared to earlier claims. Then they add to their sample six clusters 
from the literature which have strong plus weak lensing analyses available. Considering 
the 10 clusters, they claim a 7a excess of the concentration parameter compared with 
the ylCDM pred i ctions . 

lOkabe et al conducted a weak lensing analysis of a sample of 30 X-ray lu- 

minous galaxy clusters and found a mean concentration Cvir = 3.48^i['^5 for clusters 
with Mvir ~ 10^^ Mq, displaying no over-concentration. However, it is worth not- 
ing that the inclusion of strong lensing data is important for an unbiased value of the 
concentration parameter. Indeed, background galaxy catalogues may have an uncon- 
trolled degree of contamination from unlensed cluster members. This will dilute the 
she ar profile in the inne r region and bias low the inferred concentration parameter. 

ISereno et al.l 1 2010al ) investigated the over-concentration problem on a sample of 10 



strong lensing clusters. They derive which elongation along the line of sight is needed 
by the data in order to be compatible with the mass-concentration relation. Half of the 
clusters of their sample support the expectation from ylCDM simulations, being triaxial 
haloes with a strong orientation bias. The other half would fit the mass-concentration 
relation only if they were characterised by a filamentary structure extremely elongated 
along the line of sight, which the authors consider very unlikely considering standard 
scenarios of structure formations. 
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lOeuri et al ] (I2OI2I ') studied the mass distribution of 28 galaxy clusters using strong 
and weak lensing observations. The sample is made of clusters selected via their 
strong lensing properties: the Sloan Giant Arc Survey. They found the inferred mass- 
concentration relation for these clusters to be in reasonable agreement with the sim- 
ulations for very massive haloes (Mvir ~ 10^^h~^ M0). However, they found that the 
observed concentrations are much higher than theoretical expectations for less massive 
haloes (Mvir ~ lO^/i"^ M©), even after taking into account the mass dependence of 
the lensing bias. 

We have seen that, even after accounting for pro jection biases, an o ver-concentration 
problem remains in some lensing clusters (see also ICralla et all l201lh . 

Another problem concerns the number of over-concentrated clusters one can find 
in a given survey. This is an open question worth investigating in detail. We have seen 
that Abell 1689 by itself does not pose a severe challenge to the ylCDM model, but 
such high values of the concentration parameter appear to be common in the combined 
strong plus weak lensing analyses of massive clusters. 

We now come back to the Einstein radius, which provides a relatively model- 
independent measure of the mass density of a cluster core, and is observationally 
easier to infer than the concentration parameter. To date, the largest Einstein ra- 
dius has been observed in MACS J0717.5-I-3745, hereafter MACS J0717 jZitrin et all 
I2OO9I : iLimousin et aLlEoij ). Its effective Einstein radius was estimated to be 55±3" 
for a source redshift of z ^ 2.5. IZitrin et alJ (200^) claimed that the probability to find 
such a system in a /IC DM Universe is very unlikely, of the order of 10~^. Recently, 
IWaizmann et al. modeled the distribution function of the single largest mea- 

sured value of the Einstein radius in a given cosmological volume based on a Monte 
Carlo approach. Results are fitted with the general extreme value distribution. They 
showed that the large Einstein radius in MACS J0717 does not exhibit tension with 
ylCDM, even if they neglect the impact of dynamical merging which is clearly estab- 
lished for MACS J0717, an d which would allow large Einstein radii more likely to be 
found l|Redlich et al.ll2012l l. This finding is at odd s with the claims by IZitrin et al.l 
( 20091 ). and we refer the reader to Waizmann et al.l (|2012l ) for a discussion of the dif- 
ference of calculations between each studies. They concluded that, for an observed 
Einstein radius to challenge ylCDM, one should observe Einstein radii larger than 
100". 



Besides, IWaizmann et al.l (|2012l ) investigated the influence of triaxiality on the 
resulting extreme value distributions of the largest Einstein r adii, and find that it 
is ver y sensitive to very elongated objects. In agreement with lOguri and BlandfordI 
( 20091 '). they confirm that the single largest Einstein radius has not necessarily its origin 
in the most massive clusters. Instead, triaxiality, together with the haloe orientation, 
has a stronger impact than the mass of the cluster itself. 



6.5 Triaxiality & Self Interacting DM 

The asphericity of galaxy cluster sca le DM haloes demons trates that the self- interaction 
of DM particles cannot be too large. Ivbshida et al.l ( 200d ) investigated how the internal 
structure of dark haloes is affected if cold dark matter particles are assumed to have 
a large cross section for elastic collisions. It results in a cluster that is more nearly 
spherical at all radii compar ed to the coUisionless case. More recently, IPeter et al.l 



spn 

m. 



201J) and lRocha et al.l l|2012h presented cosmological simulations with self-interacting 
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dark matter (SIDM). They showed that the currently observed asphericity allow dark 
matter self-interaction cross section at least as large as ajm — O.lcm^/g. 



6.6 Dynamical Information 

A spectroscopic campaign targeting cluster members is highly relevant when studying 
a galaxy cluster and its 3D geometry. Indeed, it can give us some clues about the 
dynamical state of the system. This is particularly relevant since the more relaxed 
and unimodal the cluster, the more legitimate to describe it with a single mass clump 
(whatever the geometry considered). 

For example, in t he galajcy cluster CI 024, an apparently relaxed cluster in X-ray 
and lensing data sets, ICzoske et al.l (|200lh discovered that this cluster is composed of 
two mass clumps aligned along the l ine of sight and argued that a high-speed collision 
had taken place IjCzoske et al 1 120021 ). 



7 Conclusions: Galaxy Clusters are not Triaxial 

We have discussed that a triaxial model is clearly an important step forward to describe 
galaxy clusters more realistically than using a spherical model. 

However, as any model, it has its own limitations we aim to mention. We have seen 
in Section r6.3l that. when looking at details, haloes exhibit departures from a simple tri- 
axial geometry. Both observations and simulations show the presence of substructures 
which are not accounted for by a triaxial model for the galaxy cluster. If the substruc- 
tures are small compared to the main cluster's haloe, then the triaxial approximation 
may be accurate enough. On the other hand, unrelaxed haloes often have shapes that 
are not adequately described by ellipsoids, making shape parameters ill-defined. Indeed, 
if there is no clear dominant haloe but a superposition of sub-haloes with compara- 
ble masses, the triaxial approximation may be questionable. For example, the high 
redshift MACS clusters (Ebcling ct al. 200 7D which are found to be highly disturbed 



dEbelmg et al.i 12 004: Kar talteoe et al.l 12008: ISmith et al.ll2009l: IZitrin and BroadhurstI 
I2OO9I : IZitrin et al..,200a . l2010l : iMann and Ebelinell201ll : iLimousin et al.ll2012l') 



likely to be well described by a triaxial haloe, as proposed by Serene and ZitrinI ( 20121 ). 

We note, however, that the latter study constitutes a step forward with respect to a 
spherical analysis. Therefore, we emphasise the fact that, in order to apply a triaxial 
mass model, one should concentrate on unimodal galaxy clusters with little substruc- 
tures and as close as possible to virialization. 

We have presented in this review results obtained in a triaxial framework for four 
strong lensing clusters. To draw serious conclusions regarding the clusters' properties 
and how they might evolve with time, or with any characteristic of the cluster, we 
need to study within a triaxial framework a sample of clusters as large as possible. 
Current cluster samples, which already have complementary data sets, can allow such 
an ambitious pro ject to be started: the Local Cluster Substructure Surve y (LoCuSS 
ISmith et al.ll20ld rl: the MAssive Cluster Survey (MACS libeling et al.ll200ll ): the Clus- 



^ http: / /www. sr. bham.ac.uk/locuss/ 
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ter Lensing And Supernova survey with Hubble (CLASH Postmai^^^aJ 



2012 



and 



the Dark energy American French Team sample (DAFT iGuennou et al.ll2O10l . l2Ol 

As mentioned earlier, a triaxial model is still a simplification of what a galaxy 
cluster might be. It is not clear how departures from a strict tria:xial geometry can 
bias results based on the assumption of a triaxial model. Actually, this ought to be 
tested using numericall y simulated cluste rs: mock lensing. X-ray and SZ observations 
can be generated (e.g. iRasia et al.]|2012l ) for galaxy clusters presenting different dy- 
namical states, and analysed in the same way as observational data. Comparing the 
inferred triaxial parameters to those measured on the simulated haloes is essential to 
test the methods and to quantify how the presence of substructures and deviations 
from equilibrium can bias triaxial reconstructions. This pioneering work will also pave 
the road to more refined mass models beyond the triaxial framework. 
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